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We present the Standard Model evaluation of the indirect CP violation parameter ek using 
inputs determined from lattice QCD together with experiment: |lAs|, \Vcb\, Co, and Bk- We use the 
Wolfenstein parametrization (IPctl, A, p, fj) for the CKM matrix elements. For the central value, 
we take the angle-only fit of the UTfit collaboration, and use llAsI from the Kts and Kfj _2 decays 
as an independent input to fix A. For the error estimate, we use results of the global unitarity 
triangle fits from the CKMfitter and UTfit collaborations. We find that the Standard Model (SM) 
prediction of ek with exclusive Vet (lattice QCD results) is lower than the experimental value by 
3.6(2)(j. However, with inclusive Vab (results of the heavy quark expansion), the tension between 
the SM prediction of ek and its experimental value disappears. 


CP violation in nature was first discovered in 1964 [T] . 
The neutral kaon system has two kinds of CP violation: 
one is the indirect CP violation due to CP-asymmetric 
impurity in the kaon eigenstates in nature, and the other 
is the direct CP violation due to the CP violating nature 
of the weak interaction HE]. Here, we focus on the 
indirect CP violation, which is parametrized by £k- 

^ -AjKL -» 7r7r(J = 0)) 

^ A{Ks -T 7r7r(/ = 0)) ’ 

where Ki^ and Kg are the neutral kaon eigenstates in 
nature, and / = 0 is the isospin of the final two-pion 
state. In experiment [3], 

Ek = (2.228 ± 0.011) x 10"^ x , 

4 = 43.52 ± 0.05° . (2) 


In the Standard Model (SM), the CP violation comes 
solely from a single phase in the CKM matrix elements 
[ 5 . The mixing of neutral kaons is allowed through the 
box diagrams which describe the mass splitting AMk 
and Ek EE]. In the SM, the master formula for ek is 


Ek sin 0 Asd + ^0 + Cld) 

+ O(we') + o(eor2/ri), 

where and Asd are defined as follows. 

^ ^ GlF^rriKoM^ 

® Qy/2TT‘^AMK ’ 

Asd = f]\^\Vcb? \Vcb\‘^il - p)r]ttSo{xt)il -b r) 
+ ( 1 - ^ ) {VetSoixe, xt) - r]ccSo{xc)} 


(3) 


(4) 


(5) 


where 5 ' 0 ’s are the Inami-Lim functions [8], and xi = 
mf/M^ with i = c,t. Here, r = {t]ccSo{xc) - 


2r]ctSo{xc,Xt)}/{f]ttSo{xt)}. A, p, and 77 are the Wolfen¬ 
stein parameters of the CKM matrix elements [7]. Here, 
we replace A by Vcb, using the relation |I4h| = AX^ -b 
0(A®). r]ij with i,j = c,t represents the QCD correc¬ 
tions to the box diagrams. ^0 = ImAo/ReAg represents 
the long distance effect from the absorptive part, and 
corresponds to the long distance effect from the dispersive 
part [9]. The correction terms 0{uje') and C’(Cor 2 /ri) 
are of order 10“^, which we neglect in this paper. The 
master formula of Eq. ([^ is essentially the same as that 
of Ref. [TU]. Details on how to derive Eq. (|^ from the 
SM are given in our companion paper [^. 

In Eq. (|^, the major contribution to Ek comes from 
the Bk term, and a minor contribution of about —7% 
comes from the ^0 term. The remaining contribution of 
^ld is about 2 %, coming from the long distance effect 
on Ek [HI EH. In this paper, we neglect ^ld without 
affecting our conclusion. 

The Wolfenstein parameters X,p,fj and A can be ob¬ 
tained from the global unitarity triangle (UT) fit. Here, 
we use X,p,fj from the CKMfitter [T?l [Tl] and UTfit col¬ 
laborations [milB]. They are summarized in Table |!] 

The parameters Ek, Bk, and Vcb are inputs to the 
global UT fit. Hence, the X,p,f] parameters extracted 
from the global UT fit of the CKMfitter and UTfit 
groups contain unwanted dependence on ek, Bk, and 
Vcb- Therefore, in order to determine Ek self-consistently, 
we take another input set from the angle-only-fit (AOF) 
in Ref. [IT] . Here the advantage is that the AOF does not 
use Ek, Bk, and Vcb as inputs to determine the UT apex 
parameters p and fj. The AOF gives the UT apex {p, fj) 
but not A. We can take A independently from the CKM 
matrix element I4s, using the relation: |t4s| = A-bO(A^). 
Here, the and A ^2 decays are used to set I4s ID- 

In Table |IT[ we summarize the input values for Vcb- In 
the inclusive channel, they use B —>■ Xclv, and B —>■ A^y 
decays. They also use moments of outgoing lepton en¬ 
ergy, hadron masses, and photon energy and fit them 
to the theoretical expressions which come from the op- 
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TABLE I. Wolfenstein Parameters 



CKMfitter 

UTfit 

AOF 

A 

0.22535(65) /[3] 

0.22535(65) /g] 

0.2252(9) /g] 

P 

0.131” m 

0.136(18) /[4] 

0.130(27) /[I3 

V 

o.3451«:e;i^ m 

0.348(14) /[4] 

0.338(16) /[n] 


erator product expansion (OPE). They use perturbative 
expansion in the strong coupling constant as and inverse 
heavy quark mass A/mi,. For the b quark mass, there 
are two popular schemes in the heavy quark expansion: 
the kinetic scheme [nmni and the IS scheme |1]. Here, 
we use the results of the kinetic scheme to calculate Ek 
since they are more recently updated and have somewhat 
larger errors. 

For the exclusive channel for Vet, we use the results 
of the FNAL/MILC collaboration [20]. They use lat¬ 
tice QCD to calculate the form factor J-{w) of the semi- 
leptonic decay B —)■ D*iP at zero recoil {w = 1). The 
lattice measurements are done over the MILC Nf = 2-1-1 
asqtad gauge ensembles [21] . They use the Wilson clover 
action for the heavy quark with the Fermilab interpre¬ 
tation [22]. They combine the lattice results for T{w) 
with the HFAG average [23] of experimental values of 
-F(l)|i7EM||lch| to extract Vet- Here, |^em| represents 
small enhancement factors which come from electromag¬ 
netic corrections. 


TABLE 11. Inclusive and exclusive 11461 in units of 10~®. 
Here, Kin. represents the kinetic scheme in heavy quark ex¬ 
pansion, and IS the IS scheme. 


Inclusive (Kin.) 

Inclusive (IS) 

Exclusive 

42.21(78) /[H] 

41.96(45)(07) /g] 

39.04(49)(53)(19) /[20] 


For the kaon bag parameter Bk, we use the FLAG av¬ 
erage [24] and the SWME results [25] which are summa¬ 
rized in Table m FLAG combines several lattice results 
for Bx with Nf = 2-1-1 [25H25] to obtain the average. 
FLAG uses the Bk result of the SWME collaboration 
[26] . which is not much different from the most up-to- 
date value [2S] that we use in this analysis. The BMW 
calculation [23] has the smallest error, and it dominates 
the FLAG average. The SWME result [2S] has a larger 
error, and its value deviates most from the FLAG aver¬ 
age. 


TABLE III. Bk 


FLAG 

SWME 

0.7661(99) /[24] 

0.7379(47)(365) /[25] 


The RBC/UKQCD collaboration provides lattice re¬ 


sults for ImA 2 and in Ref. [30]. Here, we use their 
result of ^0 = —1.63(19)(20) x 10“"^ obtained using the 
experimental value of e'/e. 

The factor -qu is given at next-to-leading order (NLO) 
in Ref. [10]. The factor is given at next-to-next-to- 
leading order (NNLO) in Ref. [31]. The factor q^c is 
given at NNLO in Ref. [31] ■ In Ref. [33], they claim 
that the error is overestimated for the NNLO value of 
qcc given in Ref. [32]. Hence, in order to check the claim, 
we follow the renormalization group (RG) evolution for 
qcc described in Ref. [31] to produce the NNLO value of 
qcc- The results are summarized in Table [IV] In this ta¬ 
ble, note that the results are consistent with one another 
within the systematic errors, but our qcc value is essen¬ 
tially identical to that of Ref. [33] . Details of the SWME 
result are explained in Ref. [9]. In this paper, we use the 
SWME result for qcc to obtain Ek- 


TABLE IV. Results of qcc at NNLO. 


Collaboration Value 

Ref. 

Brod and Gorbahn 1.86(76) 

m 

Buras and Girrbach 1.70(21) 

m 

SWME 

1.72(27) 

m 

The input values for qij that we use in this 
summarized in Table [V] 

paper are 


TABLE V. QCD corrections 


Input 

Value 

Ref. 

^CC 

1.72(27) 

m 

qtt 

0.5765(65) 

m 


0.496(47) 

m 

The remaining input parameters are summarized in 
Table El] 

TABLE VI. Other Input Parameters 


Input 

Value 

Ref. 

Gf 

1.1663787(6) x 10"® GeV^ 

m 

Mw 

80.385(15) GeV 

m 

rndme) 

1.275(25) GeV 

m 

mt{mt) 

163.3(2.7) GeV 

m 

9 

43.52(5)° 

m 

m^o 

497.614(24) MeV 

m 

AMk 

3.484(6) X 10"^^ MeV 

m 

Fk 

156.1(8) MeV 

m 


Let us define as the theoretical evaluation of \ek\ 
obtained using the master formula Eq. (|^ directly from 
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the SM. We define as the experimental value of |e_R-| 
given in Eq. ([^ . We define Ae^ as the difference between 


and e^: 


Ac- — c-^P c-' 

iXek = Sk ~ ^ 


-K 


( 6 ) 

Here, we assume that the theoretical phase 9 in Eq. (§ 
is equal to the experimental phase 4>e in Eq. § M- 

we present results for obtained using 


In Table 


VII 


the FLAG average for Bk together with Vcb in both in¬ 
clusive and exclusive channels. The corresponding prob¬ 
ability distributions for and are presented in 
Fig.[T] for the AOF case. The corresponding results for 
Aek are presented in Table |V^ 


TABLE VII. in the unit of 1.0^ 10 
the FLAG average for Bk in Table 


III 


Here, we use 
The input methods 
of CKMfitter, UTfit, and AOF represent different inputs for 
the Wolfenstein parameters A, p, fj. 


Input Method 

Inclusive Vcb 

Exclusive Vcb 

CKMfitter 

2.17(23) 

1.62(18) 

UTfit 

2.18(22) 

1.63(18) 

AOF 

2.13(23) 

1.58(18) 


TABLE VIII. Aek- Here, we use from Table 


obtain a by combining the errors of and in quadra¬ 
ture. 


VH 


We 


Input Method 

Inclusive Vcb 

Exclusive Vcb 

CKMfitter 

0.24ct 

3.4cr 

UTfit 

0.20ct 

3.4(7 

AOF 

0.44cr 

3.6(7 


From Table VIII we observe no tension between 
and with inclusive Vcb- 


However, from Tables VH and VIII we find that £^ 


with exclusive Vcb is only 71% of . For this case, with 
the most reliable input method (AOF), Aek is 3.6 (t. The 
largest contribution in this estimate of £|^ that we have 
neglected is ^ld ~ 2%. Hence, the neglected contribu¬ 
tions cannot explain the gap Aek of 29% with exclusive 
Vcb- Hence, our final results for Aek are 


Aek = 3.6(2)ct 
Aek = 0.44(24)cr 


(exclusive Vcb) (7) 
(inclusive Vcb) (8) 


where we take the AOF result as the central value and 
the systematic error is obtained by taking the maximum 


difference among the input methods in Table VIII 


In the case of the FLAG Bk, the BMW result of Bk 
[29] dominates the FLAG average, and the gauge en¬ 
sembles used for the BMW calculation are independent 
of those used for the exclusive Vcb m- Hence, we as¬ 
sume that we may neglect the correlation between the 



1 1.5 2 " 2.5 


(a) ek comparison (exclusive V^b) 



(b) ek comparison (inclusive Vcb) 


FIG. 1. Probability distribution for efH (blue dotted line) 
and (red solid line). The x-axis (y-axis) is the ek value 
in the unit of 1.0 x 10“® (the probability distribution). Here, 
we use the FLAG Bk and the AOF parameters. 


FLAG Bk and the exclusive Vcb- However, the SWME 
calculation of Bk in Ref. [3S| shares the same MILG 
gauge ensembles with the exclusive Vcb determination 
in Ref. [20]. Hence, there exists a substantial correla¬ 
tion between the SWME Bk and the exclusive Vcb- We 
introduce -1-50% correlation and —50% anti-correlation 
between the SWME Bk and the exclusive Vcb and take 
the maximum deviation from the uncorrelated case as the 
systematic error due to the unknown correlation between 
them. Details of this analysis are explained in Ref. [3]. 
However, this analysis shows that the size of the ambi¬ 
guity due to the correlation between the SWME Bk and 
the exclusive Vcb is much larger than the systematic er¬ 
ror in Aek with the FLAG Bk- Therefore, we use the 
results obtained with the SWME Bk only to cross-check 
those obtained with the FLAG Bk [H]- 

In Fig. we present Aek /o’ as a function of time start¬ 
ing from 2012. In 2012, the RBG/UKQGD collaboration 
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FIG. 2. Recent history of Aex along with theoretical 
progress. 
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reported in Ref. [30]. In addition to this, using the 
LLV average for Bk |35|, we reported Aejy = 2.7(2)cr 
in Ref. |3Sj in 2012. In 2014 FNAL/MILC reported an 
updated Vcb from the exclusive channel |5D|. Using the 
FLAG average for Bk |21| and the NNLO value of rict 
pTT] . we reported the updated Asjy = 3.3(2 )ct in Ref. |37|. 
In Ref. |9], we investigate issues in the NNLO calculation 
of rycc ISH 132] , and in this paper we use the SWME re¬ 
sult in Table IV to report the updated Ae^y = 3.6(2)cr in 
Eq. 0. 

In summary, we find that there is a substantial 3.6(2)cr 
tension in ek between the experiment and the SM the¬ 
ory with lattice QCD inputs. For this claim, we choose 
the angle-only fit (AOF), the exclusive Vcb (lattice QCD 
results), and the FLAG Bk (lattice QCD results) to de¬ 
termine the central value. The systematic uncertainty 
is obtained by taking the maximum deviation from the 
central value by choosing other input methods from the 
global fits of the CKMfitter and UTfit. We choose the 
AOF method to determine the central value because the 
Wolfenstein parameters of AOF do not have unwanted 
correlation with ek, Bk, and lUchl- However, the ten¬ 
sion disappears in the case of inclusive Vcb (results of the 
heavy quark expansion based on the OPE). 

we present the error budget for for 


In Table 


IX 


the central value. This is obtained using the error prop¬ 
agation method explained in Ref. |5]. From this error 
budget, we observe that Vcb dominates the error in e|U. 
Hence, it is essential to reduce the error of Vcb as much 
as possible (see also Refs. [3H1I33])- To achieve this goal, 
there have been a lot of on-going efforts in the lattice 
community 


It is true that there is an issue with the convergence 
of the perturbative expansion of rjcc [32] . This could be 
resolved with lattice QCD calculations such as those en¬ 
visioned by the RBC/UKQCD collaboration [Ilj . 

We expect that our results for ek would be consistent 


TABLE IX. Fractional error budget for obtained using 
the AOF method, the exclusive Vcb, and the FLAG Bk- 


source 

error (%) 

memo 

Vcb 

40.7 

FNAL/MILC 

V 

21.0 

AOF 

Tjct 

17.2 

c — t Box 


7.3 

c — c Box 

P 

4.7 

AOF 

mt 

2.5 


Co 

2.2 

RBC/UKQCD 

Bk 

1.6 

FLAG 

me 

1.0 



with those from a global UT analysis, such as that in 

Ref. [35]. 
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